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Abstract 

Background: /Kuto\mmur\e thyroid disease (AITD) comprises diseases including Hashimoto's thyroiditis and Graves' disease, 
both characterized by reactivity to autoantigens causing, respectively, inflammatory destruction and autoimmune 
stimulation of the thyroid-stimulating hormone receptor. AITD is the most common thyroid disease and the leading form of 
autoimmune disease in women. Cytokines are key regulators of the immune and inflammatory responses; therefore, genetic 
variants at cytokine-encoding genes are potential risk factors for AITD. 

/Wef/joofs; Polymorphisms in the IL6-M4 QIC (rsl 800795), TNFA-i09, G/A (rsl 800629), //./B-511 CAT {rsl6944), and IFNGRl-56 
T/C (rs2234711) genes were assessed in a case-control study comprising 420 Hashimoto's thyroiditis patients, 111 Graves' 
disease patients and 735 unrelated controls from Portugal. Genetic variants were discriminated by real-time PGR using 
TaqlVlan SNP genotyping assays. 

Results: A significant association was found between the allele A in INFA-iOS G/A and Hashimoto's thyroiditis, both in the 
dominant (OR = 1.82, CI = 1.37-2.43, p-vo/ue = 4.4x10"^) and log-additive (OR = 1.64, CI = 1.28-2.10, p-vo/ue = 8.2x10"^) 
models. The allele C in //.6-174 G/C is also associated with Hashimoto's thyroiditis, however, only retained significance after 
multiple testing correction in the log-additive model (OR = 1.28, CI = 1.06-1.54, p-i/a/ue = 8.9x10"^). The group with Graves' 
disease also registered a higher frequency of the allele A in TNFA-308 G/A compared with controls both in the dominant 
(OR =1.85, CI = 1.19-2.87, p-i/o/ue = 7.0x10"^) and log-additive (OR =1.69, CI = 1.1 7-2.44, p-i/o/ue = 6.6x10"^) models. The 
risk for Hashimoto's thyroiditis and Graves' disease increases with the number of risk alleles (OR for two risk alleles is, 
respectively, 2.27 and 2.59). 

Conclusions: This study reports significant associations of genetic variants in TNFA and IL6 with the risk for AITD, 
highlighting the relevance of polymorphisms in inflammation-related genes in the etiopathogenesis of AITD. 
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Introduction 

Autoimmune diseases, which affect 5 to 7 % of tlie population 
and are frequendy responsible for severe disability, represent a 
major cause of chronic illness and an important issue in general 
healthcare [1]. Autoimmune thyroid diseases (AITD), including 
Hashimoto's thyroiditis (HT) and Graves' disease (GD), are the 
most common organ specific autoimmune disorders [2]. 



HT and GD stand for the major portion of clinical presentations 
in a wide range of thyroid autoimmune conditions, which 
culminate in thyroid dysfunction [3]. HT is a T ceU-mediated 
organ-specific autoimmune disease characterized by lymphocytic 
infiltration that leads to thyroid cells loss, clinically expressed by 
hypothyroidism [4]. In contrast, GD patients exhibit hyperthy- 
roidism, which is due to excessive secretion of the thyroid hormone 
induced by specific autoantibodies to the thyrotropin receptor 
(TSHR) produced by TSHR-reactive B ceUs [5]. 
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The etiology of HT and GD involves common as well as unique 
pathways. Both diseases carry thyroid reactive T cells that escape 
the tolerance process and infiltrate the thyroid. Nevertheless, 
distinct pathways lead thyroid-reactive T cells to either cause the 
death of thyroid cells in HT or their stimulation in GD [6] . Thus, 
there is presently awareness that the genetic susceptibility to HT 
and GD involves both shared and unique genes [7]. 

Several factors, including genetic, hormonal, environmental, 
and nutritional elements are involved in the initiation and/or 
development of AITD, however the pathophysiologic changes 
seen in AITD are mediated by inflammatory cytokines [3] . IL 1 P, 
IL6 and TNFot cytokines, and IFNGRl have proinflammatory 
activity and represent important facilitators of the immunologic 
process involved in HT and GD as components of the 
autoimmune response [8-14]. Various cytokines and respective 
genetic polymorphisms have been reported to be associated with 
AITD, yet the results are inconsistent [8-12,15]. 

In the present work, a case-control study was performed to 
assess if the risk to develop HT and GD is associated with the 
promoter single nucleotide polymorphisms (SNPs) ILlB-5 \ 1 C/T 
(rs 16944), TNFA-308 G/A (rs 1800629), /L6-174 G/C 
(rsl800795), and IFNGRl-56 T/C (rs2234711). This study 
reports significant associations of promoter SNPs in TNFA and 
IL6 with the risk for AITD and the combined effect of risk variants 
for HT and GD development. 

Materials and Methods 

Study population 

A total of 1266 Portuguese subjects were included in the study. 
The patients group comprised 531 individuals (mean age 
45. 4± 16.0 years). 111 diagnosed with GD and 420 diagnosed 
with HT. Patients were enrolled in the study from 2007 until 2013 
and the number obtained reflects the relative incidence reported 
for the two diseases (1GD/3.8HT) [16]. Patients' clinicopatholog- 
ical details are described in Table 1. The diagnosis of GD was 
established on the basis of clinical findings, decreased serum 
thyroid stimulating hormone (TSH) (<0.35 lU/mL), elevated 
serum free thyroxine (FT4) (> 1 .48 ng/ dL) and/ or free triidothyr- 
onine (FT3) (>3.71 pg/mL), positive serum antibodies to TSH- 
receptor (TRAb) (>1.8 lU/L), and typical ultrasound signs 
(hypoechogenicity and high perfusion). Diagnosis of HT was 
obtained based on clinical findings, positive serum antibodies to 
thyroid peroxidase (TPOAb) and/or thyroglobuhn (TgAb) (ac- 
cording to the method applied before or after March 2009 in the 
Department of Clinical Pathology of Hospital S. Joao), and 
characteristic ultrasound signs (hypoechogenicity and non-homo- 
geneous texture). 

The control group included 735 unmatched samples obtained 
from unrelated healthy blood donors (mean age 49. 2± 16.8 years). 
This group consisted of permanent residents in the catchment area 
of Hospital of S. Joao (Porto, Portugal), selected during the 
assembling of the EpiPorto cohort [17]. Enrollment of participants 
was performed under approval of Hospital of S. Joao ethic 
committee and included written informed consent for data and 
DNA usage. 

SNP genotyping 

Patients and controls genomic DNA was isolated from blood 
using standard proteinase K digestion with phenol/chloroform 
extraction. SNPs ILlB-511 C/T (rsl6944), TNFA-308 G/A 
(rsl800629), /L6-174 G/C (rsl800795), and IFNGRl-56 T/C 
(rs2234711) were genotyped using TaqMan Pre-Designed SNP 
Genotyping Assays (Applied Biosystems, Carlsbad, USA). PGR 
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amplification and allelic discrimination were performed according 
to product specifications with the ABI 7500 Fast real-time PGR 
system (Applied Biosystems, Carlsbad, USA). Optimization of 
each TaqMan Assay was performed using controls of known 
genot\p(' selected through DNA sequencing. Cases and controls 
were randomized during genotyping and 5% of the samples were 
genotyped in duplicate to assess the genotyping error rate 
(genotype concordance was 100%). 

Statistical analysis 

Genotype frequencies of all SNPs were obtained using SPSS 2 1 
(IBM SPSS Statistics). Compliance of alleles at individual loci with 
the Hardy-Weinberg equilibrium was measured at the level of the 

control population using a test (level of significance set to p- 
value <0.05) implemented in the SNPassoc 1.6-0 package in R 
[18]. 

Comparison of genotype frequencies between groups defined by 
status (HT and GD patients versus controls) was assessed by 
unconditional logistic regression (level of significance set to p-value 
<0.05) with the SNPassoc package implemented in R [18] and 
included gender and age. Odds ratios (OR) with respective 
confidence intervals (95% CI) were calculated considering the 
genotypic, dominant and log-additive (per allele) models of 
inheritance. The adjustment for multiple testing was performed 
by the false discovery rate (FDR) method [19]. 

Results 

In the control group, the frequencies of all SNPs did not deviate 
significantiy from those expected under Hardy-Weinberg equilib- 
rium ip-value for ILl 8-511, TNFA-30S, /L6-I74 and IFNGRl- 
56 is, respectively, 0.059, 0.122, 0.491 and 0.347). The genotyping 
success rate of all SNPs was lOO'lo in the- control and GD groups 
and >99% in the HT group. The frequency of genotypes and 
association parameters are summarized in Table 2. 

The TNFA-TiOi A allele is significantiy associated with HT risk 
in the log-additive (OR =1.64, CI = 1.28-2.10, ^-w/- 
M« = 8.2xl0"') and dominant models (OR=1.82, CI=1.37- 
2.43, p-value ^4Ax\0~% The frequency of tiie TNFA-308 A 
allele was also higher in GD than in controls in both log-additive 
and dominant models (OR = 1.69, CI = 1.17-2.44, p-val- 
Me = 6.6xl0"^ and OR=1.85, CI= 1.19-2.87, p-val- 
M« = 7.0xlO~'\ respectively). 

A significant association was found between the ILS-ll'^ allele 
C and HT risk in tiie log-additive model (OR= 1.28, CI= 1.06- 
1.54, j!)-va/Me = 8.9x10"^) but did not reach significance for the 
dominant model (OR=1.31, CI= 1.01-1.71, jf^-wato = 0.0406) 
after FDR correction (p-value <0.0106). 

Statistically significant differences were not observed in the 
frequencies of genotypes and alleles of the ILlB-511 C/T, IL6- 
174 C/G and IFNGRl -56 T/C polymorphisms between GD 
patients and controls, and the ILlB-5\\ C/T and IFNGRl-56 
C/T polymorphisms between HT patients and controls. 

The joint effect of genetic polymorphisms was assessed 
considering the high-risk genotypes TNFA A carrier and IL6 C 
carrier (Table 3). The reference group contains no high-risk 
genotypes, group 1 contains individuals with one of the high-risk 
genotypes and group 2 contains individuals with two of the high- 
risk genotypes. Individuals carrying two high-risk genotypes 
present an increased risk for GD [OR (CI) = 2.59 (1.46-4.59), p- 
value=\.2x\0-'^\ and for HT [OR (CI) = 2.27 (1.53-3.38), 
t;a/Me = 5.4x10"^. 



Discussion 

Our results report significant associations of genetic variants in 
TNFA and IL6 witii HT, and in TNFA witii GD. The risk odds 
were increased whenever an individual presented more than one 
high-risk genotype. Genetic variants in proinflammatory' cytokine 
genes can have the potential to alter the regulation of the 
transcript production or function. Considering the role of 
proinflammatory cytokines in the pathogenesis of the autoimmune 
response involved in GD and HT, a change in the function or the 
quantity of a particular cytokine may lead to the initiation or 
perpetuation of the inflammatory process. The importance of the 
ILlfi, TNFot, 1L6, and IFNGRl in the immunological response is 
well established in the Kterature. Thus, SNPs in these cytokines 
and receptor may have a role in the susceptibility to HT and GD 
in which the immune response is a major feature. 

TNFa plays an important role in the initiation of an adaptive 
immune response [3]. TNFa is produced by monocytes, T cells, 
natural killer cells, and mast cells. It is involved in the up- 
regulation of the HLA class I, activation of phagocytes, induction 
of ILl, IL6, and TNFa itself and, synergistically with IFNy, 
enhancement of the HLA class II expression [11,20—22]. TNFa 
mRNA is found in thyroid tissue from patients with GD and HT 
in higher levels than in healthy subjects [23]. Our study reports a 
statistically significant association between the 7^^4-308 allele A 
both with HT and GD. This is in agreement with previous studies 
for GD [11,14,15,24-28] but not for HT [10,25,28]. It has been 
observed that the TA'ivl-308 A allele is associated with a higher 
level of the TNFA transcript, justified by the greater potency of the 
promoter region to activate the transcription [29,30], which 
suggests that the TNFA-308 A allele may have a role in the 
pathogenesis of the AITD. Nevertheless, conflicting data were 
reported concerning the genotypic and aUelic frequencies of the 
TNFA-308 G/A polymorphism in GD and HT 
[9,10,25,28,31,32]. In fact, a meta-analysis of TNFA polymor- 
phisms in GD reported an association of the TNFA-308 SNP in 
Caucasians but not in Asians [33]. The TNFA gene is located in 
the HLA class III domain of the major histocompatibility complex 
(MHC) [34] . There are studies showing tlu- genetic contribution of 
the HLA regions for AITD susceptibility [35-38]. The HLA are 
regions of strong linkage disequilibrium, thus it cannot be ruled 
out that the association of TNFA variants with HT and GD might 
not be due to polymorphisms within TNFA itself, but rather to 
variation in a linked gene [25,28]. 

Polymorphisms in the IL6 promoter region have been 
implicated in susceptibility to carotid atherosclerosis [39], multiple 
myeloma [40], and juvenile chronic arthritis [41]. IL6 is mainly 
produced by mononuclear phagocytes, under stimulation of ILl, 
TNF or Hpopolysaccharide. IL6, a Th2 cytokine, has been 
connected with the modulation of thyroid ceUs function, and its 
expression in thyrocytes correlated positively with the degree of 
lymphocyte infiltration in HT [42,43]. Furthermore, 1L6 plays a 
major role in B ceU differentiation and T cell proliferation, and the 
deregulated production of IL6 and its receptor was related with 
the pathogenesis of autoimmune diseases by inhibition of 
autoreactive T cell apoptosis [44]. Our study reports an 
association of the C allele in /L6-174 with the risk to develop 
HT but not GD. This result is consistent with the lack of 
association with GD in a Pohsh [8] and Taiwanese populations 
[1 1]. In contrast, two other studies in an Iranian [15] and Turkish 
[27] populations reported association of the IL6-\74 SNP with 
GD, albeit in the latter the association was on the threshold of 
significance. To the best of our knowledge, there are no reports 
implying the IL6-174 G/C polymorphism in HT susceptibility. 
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Table 3. Genetic risk profile for Graves' disease and Hashimoto's thyroiditis including the susceptibility variants at TNFA and ILS". 



Graves' disease Hashimoto's thyroiditis 

Nr of risl< genotypes* OR (95% CI) p-value OR (95% CI) p-value 

0 LOO*" l.OO'' 



1 1.09(0.67-1.79) 0.7282 1.48(1.10-2.00) 0.0106 

2 2.59(1.46-4.59) 1.2x10^ 2.27(1.53-3.38) 5.4x10 



Note: The number of cases and controls genotyped for each SNP differs according to their genotyping success rates. 

^The reference group (0) contains no high-risl< genotypes, group 1 contains individuals with one of the high-risk genotypes and group 2 contains individuals with two of 
the high-risk genotypes. 

^Reference estimate. Values in bold are statistically significant with a p-value <0.05. 
doi:l 0.1 371 /journal.pone.Ol 05492.t003 



Current data suggest that the C allele at the —174 position is 
responsible for lower expression of IL6, leading to lower IL6 serum 
levels [45,46]. Other studies have demonstrated that patients with 
HT have higher serum levels of IL6, inversely correlated with the 
thyroid function [45,47,48], Nevertheless, it has been reported that 
the expression of IL6 occurs in a greater proportion of GD patients 
than in HT patients, which may correlate with the different 
lymphoid aggregate in the two AITD [43], There is no evidence 
that the /L6-174 G/C polymorphism is a causal variant for HT, 
but our results may contribute to clarify some differences in the 
pathogenesis of HT and GD. 

Polymorphisms in the ILl gene have been associated with other 
autoimmune diseases such as rheumatoid arthritis [49], inflam- 
matory bowel disease [50], and systemic lupus erythematosus [51], 
ILl [3 has pleiotropic effects, can alter cytokine production, cell 
signaling and migration, [52], Several studies have reported 
conflicting results regarding the possible role of /L7fi-511 C/T 
polymorphism in AITD. The TT genotype has been described as 
a protective variant for GD [12,53,54], however, other studies did 
not report any association between this polymorphism and GD or 
HT [10,55,56], In the present study, we also did not observe a 
significant association between the ILlB-51 1 C/T polymorphism 
and GD or HT. In agreement with our results, a recent meta- 
analysis has reported that the ILlB-511 SNP is associated with 
GD in Asians but not in Caucasians [57], 

IFNGRl encodes a class II cytokine receptor which is 
ubiquitously expressed in nucleated cells and is the receptor for 
IFNy, one of the most important Thl -related cytokines [58], 
Several functions have been attributed to IFNy, such as the 
enhancement of the expression of the HLA class I, class II and 
some adhesion molecules on thyrocytes, including intercellular 
adhesion molecule 1 (ICAMl) and lymphocyte function-associated 
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